Observation of Reduced Three-Body Recombination in a Fermionized ID Bose Gas 
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We investigate correlation properties of a one-dimensional interacting Bose gas by loading a 
magnetically trapped 87 Rb Bose-Einstein condensate into a deep two-dimensional optical lattice. 
We measure the three-body recombination rate for both the BEC in the magnetic trap and the 
BEC loaded into the optical lattice. The recombination rate coefficient is a factor of seven smaller 
in the lattice, which we interpret as a reduction in the local three-body correlation function in the 
ID case. This is a signature of correlation intermediate between that of the uncorrelated phase 
coherent ID mean-field regime and the strongly correlated Tonks- Girardeau regime. 
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The majority of experiments with quantum degener- 
ate gases have been performed in the weakly interacting 
limit, on Bose-Einstein condensates (BECs) character- 
ized by long-range phase coherence and well described 
by the mean-field Gross-Pitaevskii (GP) equation 
While the success of the GP equation in accounting for 
many experimental results has been spectacular, it has 
also led to the search for physics beyond mean-field the- 
ory. As in condensed matter physics, there is now great 
interest in highly correlated systems, where mean-field 
approaches are inapplicable. Progress toward such cor- 
related systems includes the recent observation of the 
Mott-insulator transition in BECs loaded into optical lat- 
tices [2j and the use of Feshbach resonances to increase 
interactions between atoms Q • Here we present evidence 
of strong correlations in a ID degenerate Bose gas as re- 
flected in a reduction of three-body recombination. 

The role of fluctuations and correlations in Bose gases 
increases with reduced dimension. In homogenous sys- 
tems, BEC is only possible in 3D. In 2D, a Kosterlitz- 
Thouless transition occurs, and in ID there is no finite 
temperature transition 0. By contrast, BEC is possi- 
ble in 1, 2 and 3D for trapped systems [5J]. Trapped ID 
systems with <5-function repulsive interactions are par- 
ticularly interesting, in that for high density the ground 
state is a condensate, while in the low density limit the 
ground state is a highly correlated state known as a Tonks 
gas 6]. This ground state is an example of "fermioniza- 
tion," where the repulsive interactions mimic the Pauli 
exclusion principle. Indeed, the low energy excitation 
spectrum and the correlation functions are identical to 
those of non-interacting fermions, and the many-body 
wave function of the Bose gas is equal to the absolute 
value of the fermionic wave function Q • 

For such a ID Bose gas, the degree of correlation 
depends on the ratio between two energies: the repul- 
sive energy of uncorrelated atoms at a given density, 
grim, and the quantum kinetic energy needed 
to correlate particles by localizing them with respect 
to each other on the order of the mean inter-particle 
distance d, E COI = h 2 /2md 2 . Here, g is the strength 



of the ^-function interaction, m is the atomic mass, 
and niD = l/rf is the ID density. A single parameter 
7 = E nnc /AE COI — mg/2niDh 2 entirely characterizes a 
homogeneous ID gas with repulsive short range interac- 
tions. Ref. provides the exact eigenstate solutions for 
all 7. Recently, these results have been extendend to 
harmonically trapped gases, addressing, e.g., the excita- 
tion spectrum y| , the shape of the trapped gas [2( , and 
correlation properties |lf)l | . 

The many-body ground state has two limiting forms. 
In the Tonks- Girardeau (TG) regime, where 7 3> 1, the 
ground state becomes correlated in order to minimize the 
interaction energy and the bosons become impenetrable, 
behaving like fermions as described in Q. The second 
and higher-order local correlation functions gi vanish , 
meaning that no more than one particle can be found 
at a given position. On the other hand, in the mean- 
field (MF) regime when 7 <C 1, the GP equation de- 
scribes the system well. In this regime the healing length, 
hi = ft/y/mgniD, is much larger than the mean inter- 
particle distance. Note the counter-intuitive result that 
the system reaches the correlated regime for low ID densi- 
ties, contrary to the 3D case where n^Big = ( , j/2tt) 3 ^> 1 
corresponds to the correlated regime. (Here a s is the 
zero-energy 3D scattering length and 7 is the ap- 
propriate energy ratio in 3D.) 

To probe correlations we measure three-body recombi- 
nation rates (proportional to the local third-order corre- 
lation function 53) of ID gases produced in a 2D optical 
lattice. This technique was used in Ref. 0] to demon- 
strate that there is a reduction of g% in a 3D BEC by 
a factor of six compared to a thermal gas. We observe 
a further reduction of three-body recombination in a ID 
gas compared to the 3D BEC situation. Even though 
7 ~ 0.5 for our system, which is far from the TG regime, 
this is a signature that the correlations are significant due 
to the "fermionization" of the particles. 

We realize a ID gas by confining a 3D gas sufficiently 
tightly in two directions that the radial confinement en- 
ergy hu>±_ is much larger than all other relevant ener- 
gies in the system: E unc , E COI , ksT, and the axial 
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trapping energy hto z . Since a s is much smaller than 
a± = y/h/mu)_i_ (a s /aj_ ~ 0.1 in our system), the atom- 
atom interaction strength is largely determined by a«, 
with only a small correction due to confinement |l3j, LLJ] : 
a c fr = a s /(l — 1.46 a s /^/2a±). There is no excitation 
in the radial direction and by i nteg rating over the ra- 
dial coordinates, one can show 13] that the system is 
formally equivalent to a true ID gas with interaction 
strength g = Ah 2 a e g / 'ma\, so that 7 = 2a e e/{niBa 2 L ). 
The 3D density is related to the effective ID density by 
niD = l/d = 7ra5_n 3 D- 

Our approach is to load a BEC into the ground state 
of a deep 2D optical lattice so that the BEC is divided 
into an array of independent ID quantum gases, each 
tubular lattice site acting as a highly anisotropic trap. 
Our experimental apparatus has been described else- 
where 15] . We achieve BEC with up to N = 5 x 10 5 
atoms in the (F,mj?) — (1,-1) hyperfine state of 87 Rb 
(for which a s = 5.313 nm 0). A Ioffe-P ritchard trap 
confines the atoms with initial "tight" trap frequencies 
of v x = v z = 210 Hz, and v y — 8.2 Hz, giving a peak 
atomic density of up to 3 x 10 14 cm" 3 . Before apply- 
ing the optical lattice, we adiabatically lower v x and v z 
to a "weak" trap frequency of 28 Hz, resulting in peak 
densities of ~ 5 x 10 13 cm" 3 . 

We create a 2D optical lattice from two independent, 
retro-reflected ID lattices which lie in the xy-plane and 
intersect at an angle of 80 degrees. The independent 
ID lattices are detuned from each another by 5 MHz. 
All beams derive from a Tksapphire laser operating at 
A = 810.08 nm (detuned below both 55 -> 5P tran- 
sitions at 795 nm and 780 nm), and the polarizations 
of the lattice beams are in the xy-plane. Each ID 
lattice is measure d 1 171 to be 29(1) Er deep (where 
Er = h 2 /2m\ 2 ) At each lattice site the ground 

state of the radial motion is well approximated by a gaus- 
sian wavefunction with a± — 58.5(5)nm corresponding to 
an effective ui±/2ir = 33.8(6) kHz. By observing dipole 
oscillations following a sudden, brief displacement of the 
trap center, we measure the axial frequency along the 
tubes to be uj z /2tt — 55.9(6) Hz. This frequency results 
from the combined effect of the magnetic trap and the 
dipole potential of the lattice beams along the tubes. To 
load the atoms into the lattice, the laser light is increased 
over 200 ms with an approximately half gaussian shape 
(rms width 70 ms), which is adiabatic with respect to all 
vibrational excitations. We estimate that the interaction- 
free tunnelling time from one lattice site to the next for a 
29 -Er lattice is ~150 ms. Although this is shorter than 
the time of the experiment (up to 12 s), it corresponds 
to an energy much smaller than the interaction energy in 
the tubes, and should not modify the local ID correlation 
properties [lij . 

To measure the reduction of (73 due to correlations, we 
observe the corresponding reduction in the three-body 
recombination rate coefficient. The local three-body re- 
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FIG. 1: Number of atoms as a function of time in the 2D 
lattice. The solid line is a fit to the decay as described in the 
text, and the dashed line is an extrapolation of the asymptotic 
1-body loss. 

combination rate (in either ID or 3D) is proportional to 
the cube of the local density. For 87 Rb, it is known that 
two-bod y lo sses [l2l |20| . including photoassociation at 
810 nm ]2ll]. are very small jl^. Our model, therefore, 
includes only one-body and three-body processes so that 
the total number of atoms N decays according to 

^ = - Kl N- J n 3 D dV. (1) 

We account for atomic redistribution during decay 
through the evolution of the density profile. Determina- 
tion of the three-body recombination rate coefficient 
requires an accurate estimate of the density, which we 
ascertain from a measurement of the number of trapped 
atoms as a function of time, along with a determination 
of the size and shape of the atom cloud. 

Figurc^Jshows the total number of atoms as a function 
of time t in the lattice, obtained by absorption imaging 
34 ms after release from the lattice and magnetic trap. 
We calibrate our absorption measurements by compar- 
ing the observed expansion of a released condensate to 
the known number-dependent expression for the expan- 
sion [2^]. (The inferred absorption cross section agrees 
(~10%) with one calculated from the steady state Zee- 
man sub-level distribution resulting from optical pump- 
ing.) The number of atoms in the BEC fluctuates by less 
than 20% from shot to shot. We automate the experi- 
ment to produce a BEC every minute and each data point 
is typically an average of five measurements. In order 
to minimize systematic effects due to long term thermal 
drift of the trap coils, we run current in the magnetic trap 
after the imaging such that the total time that the mag- 
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netic trap is on is the same for each measurement. While 
the atoms are in the lattice, we apply a radio frequency 
shield |20j tuned 500 kHz above the minimum of the trap 
to reduce heating without significantly increasing the loss 
of atoms trapped in the lattice. 

We measure the size of the lattice-trapped cloud in 
the xy-plane by phase contrast imaging. The initial col- 
umn density distribution of the cloud is well described 
by an integrated Thomas-Fermi (TF) profile, with radii 
of R x — 13.1(5) fiva and R y = 22.5(10) /im so that the 
observed number of atoms per tube at (a;, y) is well de- 
scribed by Mube = A/" max (l - {x/R x f - (y/i? y ) 2 ) 3 / 2 , 
where Af max = 5Nq A 2 /8ir R x R y is the number of atoms 
in the central tube and A/2 is the spacing of the tubes. 
Based on the initial total number and the measured sizes 
of the cloud, we determine A/" max = 230(40). In the xz- 
plane, we measure the size of the cloud by absorption 
imaging. The initial xz density distribution is also de- 
scribed by a TF profile, of radii R x — 15(2) ^im (in 
agreement with our phase contrast icy-measurement) and 
R z = 17(2) /im. For our parameters, the atom distribu- 
tion along the tubes (along z) is not expected to deviate 
significantly from a TF profile [g| ; indeed R z agrees with 
the ID TF value calculated based on A/" max - We note that 
the peak density is ~ 1 x 10 15 cm" 3 , which would lead 
to rapid three-body loss in a 3D system. 

We observe that the cloud slowly expands in the z- 
direction over the course of the measurement: the cloud 
expands by 5 /im in 2 s while the R x and R v radii remain 
constant. This expansion is consistent with a 1 kHz/s 
rate of energy increase, which is less than the expected 
rate due to spontaneous emission with full equilibration 
between the radial and longitudinal degrees of freedom, 
but more than the expected rate with no equilibration. 
Regardless of the mechanism, this expansion reduces the 
density only modestly during the first two seconds, when 
most of the three-body decay occurs. We take this into 
account empirically in modelling the decay. 

To model the decay using Eq. ^ we assume an overall 
3D TF density profile with gaussian radial distributions 
within each tube. In addition, for simplicity of modelling 
we assume K | D to be a constant (see below) . With these 
approximations, Eq. ^ becomes: 

^=- Kl N-a(t)Kl»N\ (2) 

where a(t) = (25/896tt 4 )(A 2 / 1 a\R x R y R z (t)f . The radii 
R x and R y are kept constant at their measured values, 
and R z (t) grows linearly in time at the measured rate 
of 2.5 /tm/s. This differential equation has an analytic 
solution which gives the total number as a function of 
time, to which we fit the data of Fig. □ 

With this analysis, we determine = 1.2(7) x 10~ 30 
cm 6 s _1 and Ki = 0.16(2) s _1 . This result is relatively 
insensitive to the specific model used for atomic spatial 
redistribution during decay, and variations among realis- 
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FIG. 2: Number as a function of time in the tight magnetic 
trap (no 2D lattice). The solid line is a fit to the decay |20||. 
and the dashed line is an extrapolation of the asymptotic 1- 
body loss. 

tic models fall within the quoted uncertainties. We at- 
tribute K\ mainly to optical pumping to states other 
than the original (1,-1) state, which are untrapped in 
the combined optical, magnetic and gravitational poten- 
tial. For the depth and detuning of our lattice, the ma- 
jority of photon scattering events return the atoms to the 
original (1,-1) Zeeman sub-level and do not contribute 
to K\. The calculated loss rate for a 29 Er lattice is 
0.17 s _1 , in good agreement with our measured K\. 

To determine the reduction of three-body recombina- 
tion in ID, we must compare to K^°. A compari- 
son in the same apparatus reduces systematic uncertainty 
due to our 15% number uncertainty. (The effect of this 
systematic uncertainty is not eliminated entirely because 
the power law dependence of dN/dt on N is different 
in ID and 3D.) We therefore repeat our experiments in 
a tight magnetic trap in the absence of a lattice, sim- 
ilar to Refs. [H EL (But see 0.) For the (1,-1) 
state, we measure [23 K$ D = 8.3(20) x 10~ 30 cm 6 s _1 
(see Fig. which is in agreement with the value of 
5.8(1.9) x 10~ 30 cm 6 s -1 measured in [l^ . 

Comparing our measurements, we find that the ratio of 
the three-body decay coefficients in ID and 3D is 0.14(9). 
This represents a factor of seven reduction in over 
<? 3D , a clear signature of correlations. 

For comparison of the observed reduction in 53 with 
theory, we calculate 7 at the center of each tube. From 
the experimentally determined density distribution we 
find at t = that 7 > 0.34, with 80% of the atoms 
having 0.34 < 7 < 0.65, and median value j m = 0.45. 
We do not expect correlations to vary significantly over 
this range of 7 0] > so the assumption of using a single 
average in the model should be reasonable. 
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FIG. 3: Comparison of our results with theoretical calcula- 
tions of the third-order correlation function g$ vs. 7. The 
solid and dotted lines represent the MF and TG limits of g : >, 
from [Toll , and the dashed line is an estimate |2^| of g% = (32 ) 3 , 
based on 32 calculated in pTj |. The measured suppression fac- 
tor /K^° is indicated by the box, where the height rep- 
resents the measurement uncertainty and the width is the 0% 
to 80% range of 7 for our system. 

In Fig. 01 we compare the measured reduction in the 
three-body loss rate coefficient with theoretical estimates 
for 33 in an interacting ID gas at T = 0. There is cur- 
rently no theory for 173 in the regime intermediate to the 
MF and TG limits (7 — 1). We have plotted an approxi- 
mate c/3 in the intermediate regime, using the expression 
g 3 ~ (32) 3 llil- and the value of g 2 for a homogenous 
system from Ref. jjj) . The value of g% is expected to be 
insensitive to T for T -C Td ~ A/" ma x^z/fcB in this range 
of 7. From the measured values of Af and uj z we estimate 
the distribution of degeneracy temperatures, finding at 
t = a peak value of ~13 kHz and a median value of 
~9 kHz. While it is difficult to measure the temperature 
in our system, the measured size at t = is consistent 
with zero-temperature TF theory, and is certainly much 
less than which during the first several seconds 

does not exceed ~ 6 kHz. 

The reduction in K^ D relative to indicates that 
we are beyond the mean-field regime, and is a signature of 
the "fermionization" of the bosonic particles in ID. The 
fundamental effects of low dimensionality on the correla- 
tion properties of a quantum Bose gas are also of practical 
interest, given the interest in the physics of "atom lasers" 
loaded in waveguides. In addition, these experiments in- 
dicate that high density, strongly correlated ID systems 
can be realized without fast decay due to three-body re- 
combination. 
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